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Abstract 

The relation between the on-shell and MS mass can be expressed through scalar and 
vector part of the quark propagator. In principle these two-point functions have to 
be evaluated on-shell which is a non-trivial task at three-loop order. Instead, we 
evaluate the quark self energy in the limit of large and small external momentum 
and use conformal mapping in combination with Pade improvement in order to 
construct a numerical approximation for the relation The errors of our final 
result are conservatively estimated to be below 3%. The numerical implications of 
the results are discussed in particular in view of top and bottom quark production 
near threshold. We show that the knowledge of new 0{af) correction leads to a 
significant reduction of the theoretical uncertainty in the determination of the quark 
masses. 
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1 Introduction 



In higher order calculations there is in general an ambiguity in the prediction for the 
physical quantities which can be traced back to the adopted renormalization scheme. 
Different renormalization conditions imply different numerical values for the parameters 
of the underlying theory. Very often it is useful to convert them from one scheme into an 
other in order to compare the final predictions for the observables in both schemes. 

In Quantum Chromodynamics (QCD) practical calculations are very often performed 
in the modified minimal subtraction (MS) scheme |^ leading to the definition of the 
so-called short-distance MS mass. The MS mass occupies a distinguished place among 
various mass definitions. First, it is a truly short distance mass not suffering from non- 
perturbative ambiguities. Second, the MS mass proves to be extremely convenient in 
multi-loop calculations of mass-dependent inclusive physical observables dominated by 
short distances (for a review see |^). On the other hand the experiments often provide 
masses which are tightly connected to the on-shell definition. Thus, conversion formulae 
are needed in order to make contact between theory and experiment. The two-loop rela- 
tion between MS and the on-shell definition of the quark mass has been obtained in ^ 
and has been confirmed in p. Until recently the accuracy of this equation was enough 
for the practical applications. Meanwhile, however, new computations have become avail- 
able which require the relation between the MS and on-shell mass at 0{a^) in order to 
perform a consistent analysis. The necessity of an accurate determination of the quark 
masses, especially those of the top and bottom ones, is demonstrated by the following two 
examples. 

The main goal of the future B physics experiments is the determination of the Cabibbo- 
Kobayashi-Maskawa matrix elements which will give deeper insight into the origin of 
CP violation and possibly also provides hints to new physics. In particular the precise 
measurement of Vcb is very promising. It is determined from semileptonic B meson decay 
rates. Thus it is desirable to know the bottom quark mass as accurately as possible as it 
enters already the Born result to the fifth power. 

One of the primary goals of a future electron-positron linear collider or muon collider 
will be the precise determination of the top quark properties, especially its mass, Mt. 
In hadron colliders like the Fermilab TEVATRON or the Large Hadron Collider the top 
quarks are reconstructed from the invariant mass of the W bosons and the bottom quarks. 
On the contrary in lepton colliders it is possible to determine the top quark mass from the 
line shape of the production cross section a{e^e~ — > tt) close to the threshold. Simulation 
studies have shown that an experimental uncertainty of about 100 MeV in the top mass 
determination can be achieved [0. Thus also from the theoretical side the ambiguities 
have to be controlled with the same precision. 

Threshold phenomena are conveniently expressed in terms of the pole mass, M, of 
the corresponding quark as this parameter naturally appears in the equations which have 
to be solved. The pole mass is a gauge invariant, infrared-finite and a renormalization 



scheme independent quantity [pi P, lIOl |TT|. However, various calculations |T2[ have shown 
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that higher order corrections have an significant infiuence on the predictions for the quark 
production cross section. It has been shown that this is connected to the definition of 
the quark mass used for the parameterization of the cross section A common 



remedy is the introduction of a new mass parameter which is supposed to be insensitive 
to long distance phenomena. One thus often speaks about short-distance masses. The 
relation of the new mass parameter to the pole mass is used in order to re-parameterize 
the threshold phenomena. On the other hand a relation of the new quark mass to the 
MS mass must be established as it is commonly used for the parameterization of those 
quantities which are not related to the threshold. 

Recently several groups evaluated next-to-next-to-leading order results for the cross 
section a{e^e~ — > ti) close to threshold [1^, |15|, |1^, |1^, In order to establish a 



consistent relation to the MS mass the three-loop relation between the MS and the on- 
shell mass is needed. This will be discussed in Section |^. 

The bottom quark mass can be obtained from sum rules for the T mesons [|l^, . 



The current strategy can be schematically outlined as follows pi], E^, |2^. Moments of the 



current correlator involving bottom quarks are determined with the help of experimental 
results on the electronic decay width and the masses of the T mesons. They are compared 
with the moments computed with the help of non-relativistic QCD. Due to the strong 
dependence on the bottom quark mass a precise determination is possible. Recently a 
next-to-next-to- leading order computation has become available p^ , [25[| . Again, for a 
consistent determination of the MS mass the order relation to the on-shell mass is 
necessary. 

A somewhat different approach for the determination of both the charm and bottom 
mass has been performed in The lower states in the heavy quarkonium spectrum 
have been computed up to order which again demands for the three-loop MS-on-shell 
mass relation in order to evaluate the MS mass. 

Another application where the three-loop relation between the on-shell and the MS 
mass is needed is connected the recent evaluation of the quartic mass corrections at order 
to the cross section a{e^e~ hadrons) E7\, ESI. For conceptual reasons the calculation 



has been performed using the MS mass. If, however, the result should be expressed in 
terms of the pole mass the corresponding relation has to be known to order a^. Note the 
Born result expanded for small quark masses doesn't have a quadratic term. Thus in this 
case the mass relation is only needed up to order and the quartic corrections require 
for the first time the 0{al) terms. 

In the present article we describe the calculation of the 0{al) correction to the mass 
relation as well as the implications of the result. A short version containing the main 
results has been published in 

The outline of the paper is as follows. In the next section the notation is introduced. 
Afterwards in Section | the method we use for the calculation is described in detail. Its 
power is demonstrated in Section ^ at the one- and two-loop level where a comparison 
with the exact results is possible. The tools needed for the three-loop calculation are 
presented in Section |^. In particular we collect all difficult integrals which are needed for 
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the computation of the fermion propagator. The three-loop results for the mass relation 
are finally discussed in Section || and their implications to quark production processes are 
mentioned in Section |^. Our conclusions are presented in Section ^ 



2 Notation 

This section is devoted to set up the notation. Throughout this paper the bare (unrenor- 
malized) mass is denoted by and the MS and on-shell renormalized ones by m and M, 
respectively. Their connection is given by the following equations 

m(/i) = Zmm^ = Zm{lj)M , (1) 

where is finite and has an explicit dependence on the renormalization scale /i. The 
main goal of this paper is the computation of Zm up to order a^. Therefore three-loop 
corrections to the fermion propagator have to be considered. A convenient variable in 
this context is 

^ = TP' <2) 

where q is the external momentum of the fermion propagator. 
The inverse quark propagator is denoted by 

[S%{q)y" = ^[m'[l-TPs)-i (l + , (3) 

where the functions Eg and depend on the external momentum g, the bare mass imP 
and on the bare strong coupling constant a^. The renormalized version can be cast in the 
form 

{SF{q))-' = i[{M-i) Sv{z) + M {zMSs{z) - Sv{z))] , (4) 

withQ 

SV{Z) = Z2{l+YPy), 

Ss{z) = Z2Zm{l-TPs). (5) 

Z2 denotes the wave function renormalization in the MS scheme which is sufficient for 
our considerations. Note that the functions 5*5 and Sy are MS quantities which later on 
are expressed in terms of the on-shell mass. The two-loop relation between m and M is 
enough to do this at order a^. 

It is convenient to write the functions Ss/v in the following way: 



Ss,v = l + T.S'^v['^] ^ (6) 



's/v 

n>l ^ " 



a. 



^Note that in contrast to the quantities defined in the wave function renormalization for functions 
Ss/v is still defined in the MS scheme. 
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where the quantities -S^/y exhibit the following colour structures (the indices S and V are 
omitted in the following): 



^(3) 



= CpSpp -\- CpCaSfa ~^ CpTuiSpL ~^ CpTSpH , 
= CpSppp + CpCaSppa + CfCa^faa + CpTniSppi + CpTSppH 
+ CpCaTtiiSpal + CpCaTSpah + CpT'^nfSpLL + CpT'^riiSpLH 
+ CpT'^SpHH . (7) 

The same decomposition also holds for the function Zm- In (|^) ni represents the number of 
light (massless) quark flavours. Cp and Ca are the Casimir operators of the fundamental 
and adjoint representation. In the case of SU (Nc) they are given by Cp = {N^ — 1)/ C^Nc) 
and Ca = Nc. The trace normalization of the fundamental representation is T = 1/2. 
The subscripts F, A and L in Eq. (0) shall remind on the colour factors Cp, Ca and Tni, 
respectively. H simply stands for the colour factor T. 

At one- and two-loop order the results for Zm read [||, H 
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^^^(/^) = Y^-yC2-^C3 + 3C2ln2 + -U + ^^M ^ -0.51056 
1111 1 . 3 . 3 . . . 185 , 11 



71 1 13 1 
= ^ + 2^2 + ^U/ + 8?M ^ 1-56205, 



^^''(Z^) = ^-C2 + ^1,m + IiIm ^ -0.15535, 



with l^M = In /i^/M^ where after the approximation signs the choice /i^ = has been 
adopted. For later use we also define the propagator of the gluon. In 't Hooft-Feynman 
gauge it is given by 

DM = (9) 
+ ze 

3 Method 

A formula which allows for the computation of the MS-on-shell relation for the quark 
mass is obtained from the requirement that the inverse fermion propagator has a zero at 
the position of the on-shell mass: 

{Sp{q)r' = 0. (10) 
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At order this requires the evaluation of three-loop on-shell integrals. Currently it is 
quite unhandy to deal with such kind of diagrams as up to now the literature lacks of 
a useful description of an algorithm for their computation. Even though the technology 
is in principle available as it was needed for the computation of the anomalous magnetic 



moment of the muon 30 



We have decided to choose a different way for the practical calculation. The starting 
point is, of course, also Eq. ([To|) . Applying it to Eq. (||) leads to the condition 

f{z) = zMSs{z)-Sv{z) = for z = l. (11) 

At a given loop-order L the Eqs. are inserted and the resulting equation is solved for 
Thus Eq. (^I]) can be cast in the form 

f{z) = g{z) + zl^^ (^^y . (12) 

Our aim is the computation of g{l). Note that the individual self energies E5 and Sy 
develop infra-red singularities when they are evaluated on-shell. The proper combination 
which leads to the relation between the MS and on-shell mass is, however, free of infra-red 
problems. 

At this point a comment in connection to Eq. ([lT|) is in order. In fact also more 
involved equations could be chosen which result in more complicated expressions for g{z). 
To us the choice in (|llD appears very natural. Furthermore the one- and two-loop results 
are reproduced with rather high accuracy. We also tried various other options; the final 
results, however, remained the same. 

The strategy for the computation of g{l) is as follows |31|: Expansions for small and 



large external momentum are computed for the quark self energies. After building the 

(13) 



proper combinations needed for g{z) a conformal mapping [32| is performed 



(l+o;)^' 

which maps the complex 2;-plane into the interior of the unit circle in the cu-plane. The 
relevant point 2; = 1 is mapped to u = 1. The motivation for this conformal mapping is 
based on the observation that the application of a Fade approximation relies heavily on 
analytic properties. Actually, g{z) develops a branch cut along the real 2;-axis starting 
from z = 1. This cut is mapped through Eq. (p!3D onto the unit circle of the u plane. 



Thus by applying Eq. (|T3D the radius of convergence is enlarged. In order to ensure also 
the convergence at the boundary a Fade approximation is performed. In the variable u 
it is defined through 

ao + aiu + . . . + amUj"" 
[mn\{uj) = — ; ; . (14) 



The details for the construction of the Fade approximants can be found in |3^. Spe- 



cial care has to be taken with the logarithms ln(— g^) which occur in the results of the 
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high-energy expansions as they could destroy the convergence properties of the described 



procedure. Details on their treatment can be found in ||33 



The described procedure is applied to each colour structure occurring in g(z) sepa- 
rately. If we assume that the decomposition introduced in Eqs. (|^) and also holds for 
g{z) the following equations are obtained 



9f{z) -- 

9fl{z) 

9FFiz) 

gFAA{z] 
gFAH{z 
gFLH{z 
gFFF^Z 
gFFA{z, 
gFFliz] 
gFFH^Z 



Sv,fl{z) 



= 

= S: 



gFA{z) = Ss,fa{z) - Sv,fa{z) , 
^, , ., gFH{z) = Ss,fh{z) — Sv,fh{z) , 

'S,FFiz) — Sv,FFiz) + Z^Ss,F , 

^s,faa{z) — Sv,faa{z) , gFAiiz) = Ss,fal{z) — Sv,fal{z) , 
Sv,fah{z) , gFLiiz) = Ss,fll{z) — Sv,fll{z) , 

>S,FLh{z) — Sv,FLh{z) , gFHH{z) = Ss,FHh{z) — Sv,FHH 

'S,FFFiz) — Sv,FFf{z) + Z^Ss,FF + Z^ Ss,F , 
's,FFa{z) — Sv,FFa{z) + Z^Ss,FA + Z^Ss^F -I 

's,ffl{z) — Sv,ffl{z) + z^S, ' ^^'"^ 



= Ss,Fiz) - S] 
= Ss,fl{z 

= s, 
= s, 

) = Ss,fah{z) 

I = 



— JS,FFL{Z) — Ov,FFL{Z) -f- Z^Os,FL + Z^Ss,F , 
= Ss,FFh{z) — Sv,FFh{z) + Z^Ss,FH + Z^ Ss,F ■ 



(15) 



The colour structures FF, FFF, FFA, FFL and FFH contain next to linear terms also 
products of lower order contributions. 

It was already realized in 



33 



that the Pade procedure described above shows less 
stability as soon as diagrams are involved which exhibit more than one particle threshold. 
In our case the interest is in the lowest particle cut which happens to be for = M^. 
The Pade method heavily relies on the combination of expansions in the small and large 
momentum region. The large momentum expansion, however, is essentially sensitive to 
the highest particle threshold. Thus, if this threshold numerically dominates the lower- 
lying ones it cannot be expected that the Pade approximation leads to stable results. 
In such cases a promising alternative to the above method is the one where only the 
expansion terms for — are taken into account in order to obtain a numerical value 
at = M^. This significantly reduces the calculational effort as the construction of the 
Pade approximation from low-energy moments alone is much simpler. In practice this 
approach will be applied if the Pade results involving also the high-energy data looks 
ill-behaved. 

In the present analysis diagrams with other cuts than for = are already present 
at the two-loop level (see Fig. P which allows us to test these suggestions. Also at three- 
loop order either = or g^ = 9M^ cuts appear. Cuts involving five or more fermion 
lines are first possible starting from four-loop order. Note that cuts involving an even 
number of fermions cannot occur. 

Let us summarize the strategy for the practical computation: The scalar and vector 
part of the one-, two- and three-loop fermion propagator is computed. Afterwards we 
construct with the help of the MS renormalization constants Z2 and and the two- 
loop relation between the MS and the on-shell mass the finite functions Ss{z) and Sv{z). 
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Figure 1: Sample diagrams contributing to the fermion propagator. The fermions in the 
closed loops may either also carry mass M or considered to be massless. 

The proper combination according to Eqs. ([ISD serves as a convenient starting point to 
apply the conformal mapping ( [TB| ) and subsequently the Fade approximation in order to 
obtain numerical results for the coefficients at the scale /i^ = M^. The use of the 
renormalization group equation leads to expression for different choices of /z^. 

As our approach relies on approximations some words concerning the error estimate are 
in order. Even at three-loop order seven terms both in the small- and large-momentum 
expansion could be evaluated. Thus we restrict ourselves to those Fade results which 
include at least terms of order and 1/z'^. Furthermore we require that the difference 
of the polynomial in the numerator and denominator is not too large, i.e. we select those 
Fade results which are close to the diagonal ones. From the results we compute the 
average and from the spread the error is estimated. The details will be specified below. 
In the next Section we will demonstrate that this prescription works very well at one- and 
two-loop order. 



4 Considerations at one- and two-loop level 

Let us at this point discuss the one- and two-loop results in great detail in order to get some 
feeling about the numerical quality of our procedure. In a ffist step we consider Feynman 
gauge and demonstrate afterwards how the results can be improved by considering a 
general gauge parameter. 

In Tab. |1] the results for different Fade approximations are listed, n indicates the 
number of low-energy moments involved in the analysis, i.e. n = 6 implies the inclusion 
of terms of 0{z^). The number of high-energy terms can be obtained in combination with 



7 



n 


P.A. 


F 


FF 


FA 


FL 


FH 


5 


[4/5] 


* - 1.0010 




-0.4881 




-3.3813 


n.5687 




-0.1613 


5 


[4/6] 


-1.0006 




-0.4862 


-k 


- 3.4099 


1.5666 




-0.0811 


5 


[5/4] 


-1.0008 




-0.4870 




-3.3936 


1.5676 




-0.1600 


5 


[5/5] 


-1.0006 




-0.4862 




-3.3699 


1.5663 


k 


-0.1616 


5 


[5/6] 


-1.0002 




- 0.4884 






1.5638 


-k 


-0.1534 


5 


[6/4] 


-1.0006 




-0.4862 




-3.3884 


1.5665 




-0.1381 


5 


[6/5] 






- 0.4884 


-k 


- 3.3972 




-k 


-0.1526 


6 


[4/6] 


-1.0007 


-A- 


- 0.4888 




-3.3148 


1.5668 


-k 


-0.1608 


6 


[5/5] 


-1.0006 




-0.4910 




-3.3628 


1.5665 


-k 


-0.1561 


6 


[5/6] 


-1.0002 


■k 


- 0.4866 


-k 


- 3.4084 


1.5640 


-k 


-0.1545 


6 


[5/7] 


-1.0003 




-0.4938 


-k 


- 3.3824 


1.5646 


-k 


-0.1572 


6 


[6/4] 


-1.0007 




-0.4894 


-k 


- 3.3830 


1.5667 




-0.1382 


6 


[6/5] 




■k 


- 0.4865 


Tk- 


- 3.4208 




k 


-0.1543 


6 


[6/6] 


-1.0003 


-k 


- 0.4897 






1.5646 






6 


[7/5] 


-1.0004 








-3.3777 


1.5648 


-k 


-0.1576 



Table 1: Results for the one- and two- loop coefficients of Zm- n indicates the number of 
moments included into the analysis. Fade approximations which develop a pole for |ci;| < 1 
are represented by a dash. Those where a cancellation with the numerator takes place 
(see text) are marked by a star {i^). For the gauge parameter ^ = has been adopted. 



the order of the Fade approximant {[x/y]) and is given hy x + y + 1 — n. 

Some Fade approximants develop poles inside the unit circle {\uj\ < 1). In general 
we will discard such results in the following and represent the results by a dash. In 
some cases, however, the pole coincides with a zero of the numerator up to several digits 
accuracy. These Fade approximations will be marked by a star {-k) and will be taken into 
account in constructing our results. To be precise: in addition to the Fade results without 
any poles inside the unit circle, we will use the ones where the poles are accompanied by 
zeros within a circle of radius 0.01, and the distance between the pole and the physically 
relevant point g^/M^ = 1 is larger than 0.1. 

The exact values are given in Eqs. (|^). At ffist sight good agreement is obtained in all 
five cases. A closer look shows that at order almost all results are consistently below 
the exact value. However, the maximal deviation is far below the per mille level. Similar 
observations can be made in the case of the light-fermion corrections at order a1 which 
are proportional to CpTni. Here the deviation from the exact result amount to a few per 
mille. The spread in the Cp term at two-loop order amounts to about 1%, however, the 
exact result is not covered. The deviation is small and well below 2%. For the structure 
CfCa the error is smaller and the correct result lies inside the interval. 

The situation is less pleasant for the FH part. Strong fluctuations are observed in 
Tab. It seems that there is a significant influence from the threshold at = 9M^ which 
has the consequence that no stabilization is observed if higher order terms in the large- 2; 
expansion are incorporated into the Fade analysis. Under these circumstances it is very 
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n 


P.A. 


FH - z 




4 


[1/3] 


-0.1552 


-0.1556 


4 


[2/2] 


-0.1553 


-0.1551 


4 


[3/1] 


-0.1553 


— 


4 


[4/0] 


-0.1552 


0.0114 


5 


[1/4] 


-0.1553 


-0.1552 


5 


[2/3] 


-0.1553 


-0.1553 


5 


[3/2] 


-0.1553 


-0.1553 


5 


[4/1] 


-0.1553 




5 


[5/0] 


-0.1553 


-0.3222 


6 


[1/5] 


-0.1553 


-0.1554 


6 


[2/4] 


-0.1553 


-0.1554 


6 


[3/3] 


-0.1553 


-0.1555 


6 


[4/2] 


-0.1553 


-0.1554 


6 


[5/1] 


-0.1553 


-0.1568 


6 


[6/0] 


-0.1553 


0.0061 



Table 2: Results for the structure FH where only low-energy moments have been incor- 
porated. Only in the case of "Fif — a conformal mapping has been used. 



promising, as was already mentioned above, to ignore the high-energy terms completely 
and compute the Pade approximations including only the knowledge from z — > 0. In this 
case a conformal mapping is not mandatory and the computation can also directly be 
performed in the variable z. The results can be found in Tab. ^ where the Pade results 
which develop poles inside the unit circle are again represented by a dash. Note that 
all possible Pade approximants are included — even the simple Taylor expansions which 
correspond to the Pades with degree zero in the denominator. Both the stability and 
the agreement with the exact result is very impressive in the case where no conformal 
mapping is performed. In the case with conformal mapping the Taylor expansions show 
large oscillations. However, the other results nicely group around the exact value with 
an uncertainty in the per mille level. For the choice /i^ = the threshold at = 
is very much suppressed as compared to the one at = 9M^. This is the reason that 
2; is a very good expansion parameter for z ^ and the sole inclusion of the low-energy 
terms leads to an excellent agreement with the exact value. Thus for FH we discard the 
results presented in Tab. and take the ones of Tab. |^ where the Pade approximation 
was performed in the variable z. 

At this point some comments in connection with the Pade method are in order. The 
final values for the mean and the errors actually depend on the Pade approximants in- 
cluded into the analysis. In Tab. |1| only those results are incorporated which involve z^, 
respectively, z^ terms and corrections of order 1/z^, 1/z^ or 1/z^ from the high energy 
expansion. We also demand that the difference in the degree of the polynomial in the 
numerator and denominator is less or equal to two. If one of these conditions is loosened 
the size of the error would only slightly increase and the range would then cover the 
exact value. Furthermore one should remember that only information from small and 
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n 


P.A. 


F 


FF 


FA 


FL 


FH 


7 


[6/7] 


— 




- 0.4954 


■k 


- 3.3596 


n.5643 







7 


[6/8] 


-1.0002 


-A- 


- 0.4926 




- 3.4026 


1.5638 


-k 


-0.1542 


7 


[7/6] 


-1.0003 


-A- 


- 0.4953 


-A- 


- 3.3603 


1.5641 




-0.1555 


7 


[7/7] 


-1.0002 




- 0.4925 


-k 


- 3.4088 


1.5638 


-k 


-0.1541 


7 


[7/8] 


— 




-0.4945 




-3.3722 


— 


-k 


-0.1517 


7 


[8/6] 


-1.0002 




- 0.4926 




- 3.3998 


1.5638 


-A- 


-0.1541 


7 


[8/7] 


* - 1.0003 










* 1.5644 


-k 


-0.1288 


8 


[6/8] 


-1.0002 


-A- 


- 0.4930 




-3.3073 


1.5639 




-0.1552 


8 


[7/7] 


-1.0002 










1.5638 




-0.1552 


8 


[7/8] 




■k 


- 0.4945 










-0.1551 


8 


[7/9] 


-1.0001 




-0.4991 


-k 


- 3.3593 


1.5632 




-0.1554 


8 


[8/6] 


-1.0002 


-A- 


- 0.4943 






1.5639 




-0.1552 


8 


[8/7] 


* - 1.0003 


■k 


- 0.4936 






*1.5646 


-k 


-0.1547 


8 


[8/8] 


-1.0001 




-0.4989 




- 3.3607 


1.5632 




-0.1554 


8 


[9/7] 


-1.0002 




-0.4993 


-k 


- 3.3594 


1.5633 




-0.1557 



Table 3: Results for the one- and two- loop coefficients where terms up to order z^, 
respectively, are incorporated. 



large external momenta enter into the Fade analysis. The numerical values, however, are 
extracted for = M^. They arise from combinations of different pieces where each one 
exhibits singularities at threshold. Thus, in general one can not expect that an agreement 
up to three or more digits with the exact result can be found. 

For the construction of the results in Tab. |l] the same input information was used 
which will be available at three-loop order. At one- and two-loop level, however, more 
moments can be computed. In Tab. ^ results are shown where terms up to order z^, 
respectively, 1/z^ are incorporated. Indeed, the general tendency is a slight reduction in 
the error for all colour factors. 

As compared to Tab. |l]the situation for FH becomes quite different. The inclusion of 
more low-energy moments leads to significantly more stable results. Note, however, that 
at three-loop order the input data for n = 7 and n = 8 are not available which means 
that for the diagrams of such type we follow the strategy outlined above. 

At three-loop level the moments are only available up to order z^, respectively, 1/z^. 
The results obtained with these input data (cf. Tabs. [l| and 0) read 

z^ = -1.0005(8), z^^ = -0.49(1), = -3.4(1), 

= 1.566(5), = -0.1553(2). ^ ' 

where the error is obtained by doubling the spread of the different Fade approximants. The 
comparison with the exact result results |^ { — 1, —0.51056, —3.33026, 1.56205, —0.15535} 
shows very good agreement — except for the structure FF where our error estimate is off 
by a factor of two. We will come back to that point below after considering the general-^ 
results. 
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n 


P.A. 


F 


FF 


FA 




n 


P.A. 


F 


FF 


FA 


5 


[4/5] 


' - 0.9982 


-0.4644 


-3.4939 




5 


[4/5] 






-0.5105 




5 


[4/6] 


-0.9989 


-0.4618 


-3.2830 




5 


[4/6] 


-1.0000 




-0.5040 


-3.3521 


5 


[5/4] 


-0.9986 


-0.4596 


* - 3.2273 




5 


[5/4] 


-1.0000 




-0.5105 


-3.3531 


5 


[5/5] 


-0.9990 


-0.4673 


-3.3145 




5 


[5/5] 


-1.0000 




- 0.5209 


-3.3499 


5 


[5/6] 


-0.9996 


* - 0.4617 


-3.2547 




5 


[5/6] 


-1.0000 




- 0.5100 


* - 3.3521 


5 


[6/4] 


-0.9989 


-0.4629 


-3.2604 




5 


[6/4] 


-1.0000 




-0.5055 


-3.3517 


5 


[6/5] 




* - 0.4616 


-3.2570 




5 


[6/5] 


* - 1.0000 


* 


- 0.5099 


* - 3.3521 


6 


[4/6] 


-0.9988 


* - 0.4468 


-3.2946 




6 


[4/6] 


-1.0000 




-0.5070 


' - 3.3576 


6 


[5/5] 


-0.9989 


-0.4770 


-3.3044 




6 


[5/5] 


-1.0000 




- 0.5119 




6 


[5/6] 


-0.9995 


* - 0.4574 


-3.2613 




6 


[5/6] 


-1.0000 




- 0.5093 


* - 3.3540 


6 


[5/7] 


-0.9994 


-0.4721 


* - 3.2672 




6 


[5/7] 


-1.0000 




- 0.5114 


-3.3472 


6 


[6/4] 


-0.9989 


-0.4670 


-3.2513 




6 


[6/4] 


-1.0000 




-0.5083 


-3.3497 


6 


[6/5] 




* - 0.4561 


-3.2675 




6 


[6/5] 


* - 1.0000 


* 


- 0.5092 


* - 3.3545 


6 


[6/6] 


-0.9994 


* - 0.4518 


* - 3.2837 




6 


[6/6] 


-1.0000 


* 


- 0.5112 


* - 3.3557 


6 


[7/5] 


-0.9994 


-0.4739 


* - 3.2628 




6 


[7/5] 


-1.0000 




- 0.5115 


-3.3464 



Table 4: Results for the structures F, FF and FA where for the gauge parameter the 
values ^ = — 5 (left table) and ^ = — 2 (right table) have been adopted. 



n 


P.A. 


F 


FF 


FA 




n 


P.A. 


F 


FF 


FA 


5 


[4/5] 


' - 1.0021 




-0.4344 




-3.3930 




5 


[4/5] 


' - 1.0039 


-0.3230 




-3.3834 


5 


[4/6] 


-1.0013 




- 0.4176 




- 3.4638 




5 


[4/6] 


-1.0024 




* 


- 3.4196 


5 


[5/4] 


-1.0017 




-0.4209 




-3.4105 




5 


[5/4] 


-1.0030 


-0.2620 




-3.3970 


5 


[5/5] 


-1.0013 




-0.4473 




-3.3778 




5 


[5/5] 


-1.0022 


-0.3700 




-3.3713 


5 


[5/6] 


-1.0005 




- 0.4212 




- 3.4175 




5 


[5/6] 


-1.0009 


* - 0.2110 


* 


- 3.4006 


5 


[6/4] 


-1.0013 




-0.4286 




-3.4050 




5 


[6/4] 


-1.0023 






-3.3917 


5 


[6/5] 






- 0.4206 




- 3.4195 




5 


[6/5] 




* - 0.1970 


* 


- 3.4017 


6 


[4/6] 


-1.0014 








-3.3301 




6 


[4/6] 


-1.0025 


-0.4670 




-3.3204 


6 


[5/5] 


-1.0013 




-0.4597 








6 


[5/5] 


-1.0024 


-0.3853 




-3.3643 


6 


[5/6] 


-1.0006 




- 0.4099 


* 


- 3.4389 




6 


[5/6] 


-1.0010 


* - 0.0740 


* 


- 3.4146 


6 


[5/7] 


-1.0007 




-0.4397 


* 


- 3.3990 




6 


[5/7] 


-1.0013 


* - 0.2749 




- 3.3856 


6 


[6/4] 


-1.0014 




-0.4353 


* 


- 3.3972 




6 


[6/4] 


-1.0025 


* - 0.3001 


* 


- 3.3857 


6 


[6/5] 






- 0.4024 


* 


- 3.4692 




6 


[6/5] 






* 


- 3.4305 


6 


[6/6] 


-1.0007 




- 0.2736 


* 


- 3.2798 




6 


[6/6] 


-1.0012 


* - 0.5118 






6 


[7/5] 


-1.0008 




-0.4444 




-3.3914 




6 


[7/5] 


-1.0014 


-0.3122 




-3.3803 



Table 5: Results for the structures F, FF and FA where for the gauge parameter the 
values ^ = +2 (left table) and ^ = +5 (right table) have been adopted. 



Up to now the discussion was based on Feynman gauge only, which corresponds to 
,^ = in our notation. In the remaining part of this section we allow for a general gauge 
parameter and demonstrate how the results of Eq. (|16D can be improved. At one- and 
two- loop level only the colour structures Cp, Cp and CfCa develop a ^ dependence. The 
proper combination of the self energies which contributes to the relation between the on- 
shell and MS mass evaluated at z = 1 must be independent of the specific choice of ^. On 
the other hand, the expansion terms actually depend on the QCD gauge parameter, ^. 
As we are dealing with a finite number of terms we are left with a residual ^ dependence 
in the final result. It is clear that for extreme values of ^ any predictive power of our 
procedure gets lost as the coefficients of ^ become dominant. However, the final values 
should be stable against small variations around ^ = 0. Thus it is worth to examine the 
dependence on ^. 

In Tabs. ^ and |^ the results are shown for the F, FF and FA structures where the 
values ^ = ±5 and ^ = ±2 have been adopted^]. The one-loop result is quite stable against 

^Despite the fact that for ^ > 1 the generating functional is in principle not defined we decided to 
choose this range for the gauge parameter. 
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the variation of ^ and even for ^ = ±5 the deviations to the exact result is small. Also 
the results for CaCf exhibit a quite small variation for the different values of ^. Only the 
structure FF shows larger deviations from the exact value if ^ is varied. Especially the 
average value for ^ = +5 is off by roughly 40%. However, one also has to take into account 
that in this case the procedure seems to be less stable and the errors are larger. A closer 
look into the results for the expansion terms shows that for this choice the coefficients of 
the terms are already dominated by the ^-dependent terms whereas, for instance, for 
^ = — 5 this is not the case. 

For all colour structures the best results are obtained for the choice ^ = — 2. This can 
be understood with the help of the following considerations. As is well known the QED 



fermion propagator get infra-red finite (in the leading log approximation) in a particular 
gauge, namely for ^ = — 2 in our notation (cf. Eq. (^). Thus, it is not surprising that 
for this particular choice our approximate results for z^, demonstrate better convergence 
to the known ones. It also suggests that the use of the same gauge condition (^ = —2) 
should lead to better results at least for the QED parts of Zm at three-loop order. This 
will be discussed in detail in Section |[ 

In Fig. the structures F, FF and FA of the function —g{z) are plotted for five 
different values of ^. In each figure ten different curves are plotted — two different Fade 
results for each value of ^. Whereas for 2; 7^ 1 the single curves show a quite different 
behaviour there is the tendency to converge to the same value for z = 1 which corresponds 
to the respective coefficient of Zm- For comparison also the result for FL is shown in Fig. |] 
which has no ^ dependence at all and which agrees with the exact result within 0.3%. In 
the case of F all curves are smooth at the point z = 1 which is reflected in the impressive 
stability shown in Tabs. 0, ^ and ^. On the contrary, for most values of ^ the two-loop 
structures FF and FA exhibit a kink for z = 1. In the case of FF even the sign of 
the derivative is different for z < 1 and z > 1 and a relatively strong variation around 
2; = 1 is observed. A closer look shows that the value ^ = —2 obviously represents an 
exception. For this choice the curves exhibit the smoothest behaviour around z = 1. This 
is also reflected in the numbers presented in Tab. ^ which are most stable for ^ = —2. 
Actually the structure FF has stronger singularities at threshold for ^ = than the other 
coefficients which explains the deviation from the exact result in Eq. ([T6|). Thus it is 
suggestive to choose ^ = — 2 in order to extract the results: 

z^ = -1.0000(0), z^"" = -0.51(2), = -3.35(1), (17) 

where the errors were obtained in the same way as before. We would like to stress that 
there is an impressive agreement with the exact results (cf. Eq. (|^)). 

Despite the fact that the uncertainties in Eqs. ( p!7D are quite small the total error 
would be significantly overestimated if one would naively add the contributions from the 
individual colour factors. Instead it is more promising to add the moments in a first step 
(of course, taking into account the correct colour factors) and performing the conformal 
mapping and Fade approximation afterwards. The results for the quantity 

= C'^z^^ + CpCAZ^^ + CpTmz^'^ (18) 
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Figure 2: z dependence of —g{z) for the structures F, FF and FA for different values of 
C,. The sohd curves belong to = —2. The dashed curves correspond to ^ = —5, 0, +2, +5 
(from short to long dashes). 

for different choices of ni can be found in Tabs. ^ and |^ where the gauge parameter has 
been fixed to ^ = and ^ = —2, respectively. In both cases the spread among the different 
results is very small. Using again twice the spread of the different Fade approximants as 
an estimate for the error we get for ,^ = 

/i)M9ht^Q^ = -14.4(6) , ^ -13.3(6) , 4')''*3^*(2) = -12.3(4) , 

/2),U9ht^^^ = -11.2(5) , 42),/*9/^t(4) = -10.2(3) , 4')''*^^*(5) = -9.1(3) . (19) 

Good agreement with the exact results, which are listed in the last column of the tables, 
is obtained. For ^ = —2 the following results can be extracted 

^m^ght^Q^ = -14.3(2) , 42),z.5ht(^) ^ -13.2(2) , 2^)''*3^*(2) = -12.2(2) , 
^mMght^^^ = -11.1(1) , 42),z.ght^4) ^ _io.07(3) , 4')''*3'^*(5) = -9.02(1) . (20) 



Here the error is slightly reduced and is well below 



70. 
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n 


P.A. 


n; = 


ni = 1 


ni = 2 


m = 3 


n( = 4 


ni =5 


5 


[4/5] 




-14.3431 




-13.2932 




-12.2433 




-11.1933 




-10.1432 




-9.0931 


5 


[4/6] 




-14.2673 




-13.2188 




-12.1700 




-11.1206 




-10.0703 




-9.0183 


5 


[5/4] 




-14.3816 




-13.3305 




-12.2791 




-11.2273 




-10.1750 




-9.1221 


5 


[5/5] 




-14.3188 




-13.2702 




-12.2216 




-11.1730 




-10.1244 




-9.0758 


5 


[5/6] 




- 14.7124 


-k 


- 13.5777 




- 12.4638 




- 11.3641 




- 10.2745 




-9.1925 


5 


[6/4] 


Tt 


- 14.4299 




- 13.3703 


-k 


- 12.3109 




- 11.2517 




- 10.1927 




- 9.1340 


5 


[6/5] 
















- 11.4712 


* 


- 10.3346 


* 


- 9.2264 


6 


[4/6] 




-14.2802 




-13.2344 




-12.1885 




-11.1425 




-10.0963 




-9.0499 


6 


[5/5] 




-14.3149 




-13.2662 




-12.2176 




-11.1689 




-10.1203 




-9.0717 


6 


[5/7] 




- 14.5630 


* 


- 13.4629 




- 12.3745 


* 


- 11.2942 


* 


- 10.2200 


* 


- 9.1502 


6 


[6/4] 




- 14.3676 




- 13.3163 




- 12.2647 




- 11.2130 




- 10.1610 


-k 


-9.1088 


6 


[6/5] 




-14.2190 




-13.1684 




-12.1164 




-11.0623 




-10.0043 




-8.9382 


6 


[6/6] 




- 14.2867 


* 


- 13.2389 




- 12.1910 


* 


- 11.1431 


* 


- 10.0951 


* 


- 9.0469 


6 


[7/5] 




-14.3713 




-13.3191 




-12.2668 




-11.2142 




-10.1613 




-9.1082 




exact 


-14.2287 


-13.1874 


-12.146 


-11.1046 


-10.0633 




-9.02188 



Table 6: Results for for different values of ni. The gauge parameter has been 

fixed to ^ = 0. 



n 


P.A. 


ni = 


ni = 1 


ni = 2 


ni=3 


711=4 


ni =5 


5 


[4/5] 




-14.2809 




-13.2304 




-12.1797 




-11.1284 




-10.0754 


* 


- 9.0136 


5 


[4/6] 




-14.1864 










k 


- 11.1492 


k 


- 10.0769 




-9.0194 


5 


[5/4] 




-14.3044 




-13.2494 




-12.1933 




-11.1360 




-10.0771 




-9.0164 


5 


[5/5] 




-14.2670 




-13.2184 




-12.1698 




-11.1212 




-10.0726 




-9.0245 


5 


[5/6] 


-k 


- 14.3356 


k 


- 13.2670 


k 


- 12.2014 


k 


- 11.1382 


k 


- 10.0769 




-9.0173 


5 


[6/4] 


-k 


- 14.3073 




- 13.2492 




-12.1913 




-11.1336 




-10.0760 




-9.0186 


5 


[6/5] 


-k 


- 14.3458 


* 


- 13.2720 




- 12.2035 


•k 


- 11.1387 


* 


- 10.0769 


* 


- 9.0173 


6 


[4/6] 




-14.2460 




-13.1986 




-12.1504 




-11.0983 




- 10.0940 




-9.0191 


6 


[5/5] 




-14.2625 




-13.2139 




-12.1653 




-11.1167 




-10.0677 




-9.0202 


6 


[5/6] 






* 


- 13.3253 


★ 


- 12.2248 


★ 


- 11.1461 


* 


- 10.0785 




-9.0172 


6 


[5/7] 


k 


- 14.3106 


kr 


- 13.2490 


k 


- 12.1892 


k 


- 11.1309 




-10.0737 




-9.0177 


6 


[6/4] 


k 


- 14.2928 




- 13.2394 


k 


- 12.1853 


k 


- 11.1306 




-10.0748 




-9.0179 


6 


[6/5] 










•k 


- 12.2729 


•k 


- 11.1535 


* 


- 10.0790 




-9.0175 


6 


[6/6] 




- 14.2375 


k 


- 13.1871 


k 


- 12.1322 








- 10.0819 


ir 


-9.0177 


6 


[7/5] 




-14.2909 




-13.2370 




-12.1827 




-11.1280 




-10.0731 




- 9.0177 




exact 


-14.2287 


-13.1874 


-12.146 


-11.1046 


-10.0633 




-9.02188 



Table 7: Results for for different values of ni. The gauge parameter has been 

fixed to ^ = —2. 



In summary, very good agreement with the exact results for the different coefficients 
of Zm are observed at one- and two-loop order. The results obtained from the fifth and 
sixth low-energy moment and a high-energy expansion up to 0{l/z^) suggests that also at 
order it is possible to get reliable results. The analysis for different values of ^ showed 
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Figure 3: Generic three-loop bubble diagram. In general each line may be massless or 
carry mass M. 

that the choice ^ = —2 would be the preferable one. However, also ^ = leads to very 
good results. At three-loop order the calculation for general ^ is only possible for some 
of the diagrams, namely the ones which are already present in QED. For those diagrams, 
which have stronger singularities for z = 1, ^ = —2 will be adopted. In the remaining 
cases the results are obtained for ^ = 0. In the approach where the moments are added 
in a first step our method provides for = results which are in very good agreement 
with the exact ones. 



5 Master integrals 

Before we go on to the discussion of the three-loop results the technology behind the 
practical calculation shall be discussed in more detail. The high-energy expansion can 
be reduced to the evaluation of massless propagator-type diagrams up to three loops and 
one- and two-loop vacuum graphs. The relevant tools are available since quite some time 
(see, e.g., [0). The technology needed in the limit 2; — ^ 0, however, has never been 
applied to real processes. Furthermore the discussion in the literature is not complete. In 
this Section we review the knowledge about vacuum graphs and close the gaps. 

The computation of the three-loop diagrams in the limit of vanishing external momen- 
tum reduces to the evaluation of vacuum graphs as shown in Fig. ^ with one dimensionful 
scale given by the mass, M. In principle each line may be massless or carry mass M. The 
general strategy to compute such kind of diagrams is based on the integration-by-parts 
method P^. It provides recurrence relations which are used in order to express compli- 



cated integrals in terms of simple ones and a small set of a few so-called master integrals 
for which a real integration is necessary. Those master integrals, however, only have to 
be computed once and forever. 

The recurrence relations and master integrals which are necessary for the computation 
of 0{a1) corrections to the 7 and Z boson propagator have been considered in [^. In 
this case only one difficult^ integral had to be solved. It is shown in Fig. §(e). In the 



3 



Here we mean those master integrals which cannot be solved by successive application of one- and 
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(a) (b) (c) 

Figure 5: Two- and three- loop master diagrams corresponding to the expressions in 
Eqs. (p^ ) and (|25D. The full lines carry mass M whereas the dashed ones are massless. 

technique has been extended to the case of the W boson. Altogether only three difficult 
master integrals are necessary for the evaluation of the gauge boson self energies. In 
addition to the one in Fig. ^(e) also the diagrams in Fig. ^(c) and Fig. |^(b) have to be 
taken into account. 

For the fermion propagator the variety of integrals to be computed at three-loop level 
is larger as there are more choices whether a certain line of the three-loop tadpole diagram 
shown in Fig. ^ may be massive or massless. Let us collect all master integrals which are 
needed for the computation of the moments for the fermion propagator. It is convenient 
to introduce the following function 

M,,,„,,„,(Mi,...,M6) = (21) 
two-loop formulae. 
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(M%)3D/2 (p2 + M2)-i(p2 + M|)"^(pi + M|)(p| + M|)(pi + M|)(pi + Mi) 



where the masses Mj are either equal to M or zero. D = 4 — 2£ is the space-time 
dimension. The factors of have been introduced in order to make the function M(. . .) 
dimensionless. The momenta pi are hnear combinations of the integration momenta p, k 



and /. In |^0[ the recurrence relations for those mass decompositions were derived which 



were not discussed in I^Sj. They can be used in order to arrive at the master integrals 
pictured in Fig. |[ The results can be found in |3^, |4l[| P| 



Mni(M, M, M, M, M, 0) 



2^3^.. 469 8 

h oG C4 H — 

e ^ 27 ^ 3 



CI2 



16 E 

m>n>0 



;-l)'^cos(27™/3) 



Min(M, M,0,M,0,M) 
Miii(0,0,M, M,0,M) 
Mni(M, M, 0,0, 0,0) 

Min(0,0,M, M, M, M) 
Miii(0,M, M,0,M,0) 



2(3 
2(3 

e 

2(3 
2(3 
2(3 
2(3 



- 8.2168598175087380629133983386010858249695 , 
77 

21 2 
+ 6(3 - y (4 - 4(2 In' 2 + - In^ 2 + 16Li4 ^ ^ 

+ 6(3 - 22(4 - 8(2 In' 2 + ^ In^ 2 + 32Li4 



6 


CI2 


m 










CI2 




6 











+ 6(3 - ^(4 - 4 



CI2 



(22) 



where („ is Rieman's zeta function with the values (2 = 7r'/6, (3 ~ 1.20205690 and 
(4 = 7r'^/90. Cl2(a;) is the Clausen function and Li4(x) represents the quadrilogarithm. 
The order of the expressions listed in Eq. ( p2D agrees with the one of Fig. ^(a)-(f). Note 
that the case where all six lines are massive does not occur for the three-loop QCD 
corrections to the quark propagator. 

It is interesting to note that the results corresponding to Figs. ||(a) and (e) can be 
obtained from the 0{e) part of the diagrams where the index of one of the massless lines 
is reduced to zero using the integration-by-parts technique. The resulting topology leads 
to much simpler expressions as it can essentially be written as a product of two one-loop 
integrals with an additional integration over the corresponding external momentum. For 
the other four cases the topologies presented in Fig. ^ have to be considered in order 
to reproduce the results of Eq. (^). The reason for this is that the absence of any of 
the lines immediately leads to (in general simpler) integrals which at the end result in 
different master integrals. 



"'Following standard MS practice 7^ and ln47r is discarded in the in the explicit expressions for the 
integrals. In practice this means that we multiply with a factor Q^^'^i^^)+'^iE fQj. each loop. 
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In the practical application of the recurrence relations it turns out that except the 
master integrals of Eq. (^) also the one of Fig. §(c) is needed. Furthermore the 0{D — 4) 
parts of the diagram pictured in Fig. ^(b) is required^. These pieces will be provided in 
the following. 

For their practical computation it is convenient to consider integrals which are less 
divergent. This is achieved by doubhng some of the propagators. In the case of Fig. H(b) 
we double each of the massive propagators and for Fig. ^(c) the propagators corresponding 
to line 1, 2 and 4 are raised to power two. This leads us to ultra-violet convergent integrals. 
The infra-red divergence introduced in the latter case does not spoil the evaluation. It 
manifests itself as as single 1/e pole in the final result. In a next step we introduce 
Feynman parameters and thus combine the different denominators. Then the momentum 
integrations can be carried out and one ends up with two- (Fig. §(b)), respectively, three- 
dimensional (Fig. P(c)) integral representations. In each case the leading terms of order 
1/e, respectively, e° can be used in order to check the constants which are already known. 
The expansion up to next-to-leading order in e provides the structures we are interested 
in. In both cases one ends up with a one- dimensional integration where the most difficult 
integrals look like: 



1 x" ln(x) ln(l — X + x^) 

dx 

[1 — X + X 



2\2 



n 



0,1, 



(23) 



The integration can be performed analytically and finally leads to the following results 
for the master integrals: 



Mooi(0,0,M, M,0,M) 



Mno(0,M,M,0,M, 0) 



15 1 /65 



TT 

3 



135 45 



^ + ^ + T ( - + 7tC2 ) + — + -C2 - C; 
+ 3V3CI2 



2'V 16 
763 QVStt 



32 



16 



ln^3 



35V37r3 
48 



195 15 57 
+ ^C2--C3 + -C4 

-ITT /a \ 



2 V3 



- 27V3Im 



11 



Lis 



3^3 e 



-(-lA- 2C2 + 2v^Cl2 



TT 



139 



TT 



ln^3 



17v/37r3 



72 



+ IOV3CI2 ( I ] - 6\/3Im 



Lis 



21 1 

YC2 + -C3 

-j7r/6 ' 



(24) 



These expressions were used in |1| and coincide with the results listed in a recent work E2| 



^For some diagrams it is necessary to know the 0{D — 4) part as the recurrence relations may contain 
artificial 1/{D — 4) poles. 
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As an outcome of the recurrence relations also products of one- and two-loop integrals 
may appear. This makes in necessary to evaluate also the 0{e) part of the two-loop dia- 
gram shown in Fig. ^(a). A straightforward Feynman parametrisation gives (in agreement 
with pi): 



T(M, M, M) 



(M%)« (p2 + M2)(A;2 + M2)((p + A;)2 + M^) 



3 

2^2 



9 

2e 



21 

y 



-(2 + 2^3012^2 



TT 



35V37r3 
216 



+ Cs + 6^3012 ) - eVsim 




(25) 



All other integrals which have to solved are of one- and two-loop tadpole type. The 
two- loop integrals are similar to the ones of Eq. (^), however, one or two lines are 
massless. For these kind of integrals there exists a solution valid to any order in e and for 
arbitrary exponents of the propagators in the denominator. 

A numerical check of the integrals in Eqs. (^) has been performed. A straightforward 
Feynman parameterization leads to multidimensional integrals for which a Monte Carlo 



integration routine has been used in order to obtain 5 significant digits. 

As an independent check of the results for the diagrams pictured in Fig. ^ the following 
method was used: One of the massive lines is cut in such a way that the remaining two- 
loop diagram has a threshold starting from AM"^. Then the mass of the cut line is set 
to M' and an asymptotic expansion for M' -C M is performed. The effective expansion 
parameter is (M')V(4M2). Thus is very suggestive that even for M' = M a reasonably 
fast approximation to the exact result can be expected. Actually the first 20 expansion 
terms for the diagrams in Fig. ^reproduce the first 10 digits of the exact result. Note that 
it only makes sense to apply this method if the diagram which results from the original 
one after setting the mass of the selected line to zero is computable analytically. 



6 Three-loop results 

In this section the results of order are discussed. Altogether 131 diagrams contribute. 
The evaluation in the limit z —>■ reduces to a naive Taylor expansion in the external 
momentum. Thus one ends up with three-loop vacuum integrals where the computation 
follows closely the lines of Section |^. In the limit of large external momentum the rules 
of asymptotic expansion (see, e.g., Refs. ^ for reviews) have to be applied. As a 



result for each diagram several subgraphs have to be considered which leads to roughly 
3500 diagrams to be evaluated. It is clear that this cannot be done by hand and the use 
of computer algebra becomes unavoidable. The practical calculation has been performed 
with the package GEFICOM |g6|. It calls the programs LMP respectively, EXP ^ in 



order to perform the asymptotic expansion. A recent summary of the available software 
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can be found in [^. In both limits we were able to evaluate the first seven terms in 
the expansion, i.e. terms up to 0{z^) and 0{l/z^), respectively, are available. Thereby 
Feynman gauge has been adopted. 

In principle the expansion terms of the colour structures FFF, FFA, FFL, FFH, 
FAA, FAL und FAH depend on the QCD gauge parameter, ^. The remaining three 
structures contain two closed fermion loops and thus C, drops out at an early stage of 
the calculation. A straightforward computation of the ^-dependent terms produces huge 
expressions in intermediate steps and leads quite fast to the limitations set by the soft- 
and hardware. For the abelian case there exists a simple Ward identity which allows a 



very fast calculation of the ^-dependent part of the fermion propagator 



Within dimensional regularization Eq. (^) is an exact equation which relates the ^- 
dependent part of the fermion propagator in L + 1-loop order with the full propagator at 
L-loop order. In fact, the calculation of the r.h.s. of (|26| ) for general gauge is significantly 
faster than that of the very function Sp^^ in (the simplest) Feynman gauge. This is 
because of the "factorizable" nature of the integration with respect to the momentum p. 

Using Eq. (^) we have computed the (^-dependent) QED colour structures FFF, 
FFL and FFH. To our best knowledge there is not a simple generalization of the 
identity ( pG] ) for QCD, i.e. the non-abelian contributions. Thus the results for the corre- 
sponding colour factors can only be analyzed for ^ = 0. 

Let us in a first step discuss the results for ^ = and afterwards consider the QED-like 
diagrams for general ^. 

In Tabs. || and | the results for the individual colour structures can be found where ^ = 
has been adopted for the gauge parameter. Only those results are listed which include 
at least the terms of 0{z^) and 0{l/z^). Furthermore we require that the difference of the 
degree of the polynomial in the numerator and denominator is less or equal to two. The 
Fade approximations which develop poles inside the unit plane and where no appropriate 
cancellation takes place are again represented by a dash. 

Very nice agreement between the different values in each column of Tabs. § and ^ 
can be found. There is one obvious exception in the structure FFH namely the [5/7] 
Fade approximation where the sixth moment was used for the construction. A closer look 
shows that there is a pole close to z = 1 [u = 1.105 . . ., respectively, z = 0.997 . . .) which 
explains the deviation from the other results. In order to obtain the final numbers we will 
not take this value into account. 

From Tab. |^ one observes that the structures FLH and FHH are less stable than the 
others. The reason for this is the same as for FH at two-loop order. Again the threshold 
at = is suppressed as compared to the one at = 9M^ which suggests that the 
expansion around z = should provide a reasonable expansion parameter. In Tab. TO 



the results for these two colour structures where no high energy terms were used are 
shown. Actually almost all Fade approximants performed in the variable u develop poles 
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n 


P.A. 


FFF 


FFA 


FFL 


FFH 


FAA 


5 


[4/5] 


2.0338 


-2.0665 


0.7347 


2.4876 


-16.2692 


5 


[4/6] 


2.0272 


-2.0031 


0.7341 


2.5147 


* - 16.4584 


5 


[5/4] 


2.0310 


-2.0632 


0.7371 


2.4752 


-16.3826 


5 


[5/5] 


2.0263 


-1.9981 


0.7343 


2.5057 


-16.4243 


5 


[5/6] 


2.1191 


* - 1.9318 


*0.7379 


2.5305 


* - 16.4890 


5 


[6/4] 


2.0344 


-2.0025 


0.7325 


2.5623 


-16.4381 


5 


[6/5] 


2.1091 


* - 1.8018 


0.7421 


2.5281 


* - 16.4903 


6 


[4/6] 


2.0236 


* - 2.0587 


0.7411 


2.4954 


-16.3748 


6 


[5/5] 


2.0220 


-2.0457 


0.7423 


2.4951 


-16.3949 


6 


[5/6] 


^1.9686 


* - 2.0304 


*0.7372 


2.5188 


* - 16.4554 


6 


[5/7] 


*2.0648 


* - 2.0962 


0.7601 


1.1897 




6 


[6/4] 


2.0234 


-1.9897 


0.7410 


2.4984 


-16.4185 


6 


[6/5] 


2.2334 


* - 2.0243 


*0.7392 


2.5182 


* - 16.4689 


6 


[6/6] 


^2.0548 


* - 2.0858 


0.7653 


*2.6713 


^- 16.4411 


6 


[7/5] 


^2.1493 


* - 2.1205 


0.7567 


2.3704 


-16.3435 



Table 8: Three-loop results for (X = FFF, FFA, FFL, FFH, FAA) for ^ = 0. The 
same notation as in Tab. |1] has been adopted. 



n 


P.A. 


FAL 


FAH 


FLL 


FLH 


FHH 


5 


[4/5] 




* - 1.4898 


* - 1.9890 


-0.0630 


-0.1806 


5 


[4/6] 


13.4932 


-1.3224 




-0.0708 


-0.1958 


5 


[5/4] 


^13.6711 


-1.5551 


* - 1.9764 


-0.0601 


-0.1669 


5 


[5/5] 


13.4519 


* - 1.5357 




-0.0686 


-0.1948 


5 


[5/6] 


*13.7944 


* - 1.5136 


* - 2.0122 


-0.0703 




5 


[6/4] 


13.5259 


-1.4147 




-0.0824 


-0.3390 


5 


[6/5] 


^13.6020 


* - 1.5133 


-1.9646 


-0.0690 




6 


[4/6] 


13.3826 


* - 1.5250 


* - 1.9937 


-0.0662 


-0.1939 


6 


[5/5] 


13.3533 


* - 1.5254 


* - 1.9985 


-0.0660 


-0.1930 


6 


[5/6] 


12.9396 


* - 1.5183 


-1.9767 


-0.0695 


-0.1928 


6 


[5/7] 


13.3481 




-1.9696 


-0.0486 


-0.1887 


6 


[6/4] 


13.4640 


-1.4777 


* - 1.9981 


-0.0676 


-0.2154 


6 


[6/5] 




* - 1.5173 




-0.0695 


-0.1928 


6 


[6/6] 


13.3138 




-1.9698 


-0.0319 


-0.1911 


6 


[7/5] 


13.4209 




-1.9693 


-0.0542 


-0.1880 



Table 9: Three-loop results for (X = FAL, FAH, FLL, FLH, FHH) for ^ = 0. The 
same notation as in Tab. |I| has been adopted. 



for |cij| < 1. Thus we decided not to use them for our analysis. In analogy to the two-loop 
case the results are very stable and only a small differences between the naive Taylor 
expansion and diagonal Pade results are observed. Thus in the following the results of 
Tab. |TD| are used for the structures FLH and FHH . 

It is tempting to consider also for the colour factors C|,T and CfCaT the results 
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n 


P.A. 


FFHz 


FAHz 


FLHz 


FHHz 


4 


[1/3] 


2.4135 


-1.5047 


-0.0669 


-0.1922 


4 


[2/2] 


2.4135 


* - 1.4790 


-0.0670 


-0.1922 


4 


[3/1] 


2.4065 


-1.5042 


-0.0670 


-0.1922 


4 


[4/0] 


2.3531 


-1.494 


-0.0673 


-0.1922 


5 


[1/4] 


2.4310 


-1.5085 


-0.0669 


-0.1922 


5 


[2/3] 


^2.4143 


-1.5103 


-0.0669 


-0.1922 


5 


[3/2] 


2.4364 


-1.5102 


-0.0668 


-0.1922 


5 


[4/1] 


2.4236 


-1.5080 


-0.0669 


-0.1922 


5 


[5/0] 


2.3765 


-1.4990 


-0.0671 


-0.1922 


6 


[1/5] 


2.4410 


-1.5109 


-0.0668 


-0.1922 


6 


[2/4] 


2.4473 


-1.5126 


* - 0.0669 


-0.1922 


6 


[3/3] 


2.4488 


-1.5120 


-0.0668 


-0.1922 


6 


[4/2] 


2.4474 


-1.5126 


-0.0668 


-0.1922 


6 


[5/1] 


2.4352 


-1.5105 


-0.0668 


-0.1922 


6 


[6/0] 


2.3933 


-1.5024 


-0.0669 


-0.1922 



Table 10: Results for the colour structures involving a closed heavy fermion loop. Only 
the moments for 2; — have been implemented into the Fade analysis. 



which exclusively contain the information from the expansion around z ^ Q. They are 
also displayed in Tab. |l^. There is perfect agreement with the numbers extracted from 
Tabs. P and ^ the errors, however, are significantly smaller. Thus we decided to take also 
for these colour structures the results of Tab. p!0| . 

Finally we get the following results if ^ = is adopted: 



(27) 



where the error has again been obtained by doubling the spread of the Fade approximants. 

Let us now analyze the QED-like diagrams, i.e. the structures FFF, FFL and FFH 
for arbitrary gauge parameter. In analogy to the one- and two-loop case considered in 
Section ^ we vary ^ between —5 and +5. The corresponding results can be found in 
Tabs. and |l|. 



FFF 


= 2.1(3), 


^FFA _ 


-2.0(4) , 


FFL 
m 


= 0.74(4), 


^FFH 


= 2.4(1), 


FAA _ 


-16.4(3), 


^FAL 
m 


= 13(1), 


^FAH 


= -1.51(5), 




-1.98(6), 




= -0.0669(6) 


^FHH 


= -0.1922(0), 











In the case FFH we actually prefer to use only the low-energy moments. The corre- 
sponding results for ^ = —5, —2, 0, +2 and -|-5 can be found in Tab. [1^. For all choices very 
stable results with a spread below 3% are obtained. Furthermore they are all consistent 
with each other. 

It is instructive to examine also at three-loop order the z dependence around z = 1 
for different values of ^ which is done in Fig. ^ for the ^-dependent colour structures 
FFF, FFL and FFH. The result for FLL which does not depend on ^ is plotted for 
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n 


P.A. 


FFF 


FFL 


FFH 




n 


P.A. 


FFF 


FFL 


FFH 


5 


[4/5] 


2.1093 


1.0402 


2.5059 




5 


[4/5] 


2.0599 


0.8596 


2.4949 


5 


[4/6] 


2.0926 


1.0354 


2.4652 




5 


[4/6] 


2.0447 


0.8527 


2.4978 


5 


[5/4] 


2.1082 


1.0430 


2.5013 




5 


[5/4] 


2.0580 


0.8599 


2.4865 


5 


[5/5] 


2.0923 


* 1.0443 


2.4650 




5 


[5/5] 


2.0444 


0.8549 


2.4932 


5 


[5/6] 


*1.5554 


1.0229 


3.1456 




5 


[5/6] 


2.3098 


0.8545 


2.5776 


5 


[6/4] 


2.0954 


1.0193 


2.4717 




5 


[6/4] 


2.0489 


0.8468 


2.5242 


5 


[6/5] 


2.6737 


1.0230 


2.8145 




5 


[6/5] 


2.2055 


0.8558 


2.5769 


6 


[4/6] 


2.0985 


'1.0540 


2.4565 




6 


[4/6] 


2.0434 




2.4867 


6 


[5/5] 


2.1016 


* 1.0602 


2.4248 




6 


[5/5] 


2.0369 


*0.8575 


2.4863 


6 


[5/6] 


*2.0780 


1.0016 


* 2. 4470 




6 


[5/6] 


* 2.0280 


0.8534 


2.5959 


6 


[5/T] 


*2.1475 


* 1.0204 


*2.5472 




6 


[5/7] 


* 2.0854 


0.8568 


*2.5867 


6 


[6/4] 


2.1000 


1.0283 


2.4623 




6 


[6/4] 


2.0441 


0.8537 


2.4863 


6 


[6/5] 


*2.0498 


1.0103 


*2.3488 




6 


[6/5] 


* 1.9423 


0.8535 


2.5567 


6 


[6/6] 


*2.1400 


* 1.0278 


*2.5325 




6 


[6/6] 


* 2.0777 


0.8570 


*2.5503 


6 


[7/5] 


*2.1828 


*1.0118 


*2.6241 




6 


[7/5] 


* 2.1340 


0.8569 


*3.5102 



Table 11: Results for the structures FFF, FFL and FFH where for the gauge parameter 
the values ^ = — 5 (left table) and ^ = — 2 (right table) have been adopted. 



n 


P.A. 


FFF 


FFL 


FFH 




n 


P.A. 


FFF 


FFL 


FFH 


5 


[4/5] 


1.9184 


0.6080 


2.4799 




5 


[4/5] 


0.9431 


0.4164 


2.4671 


5 


[4/6] 


1.8976 


0.6182 


2.5284 




5 


[4/6] 


1.2669 


0.4493 


2.5427 


5 


[5/4] 


1.9030 


0.6138 


2.4614 




5 


[5/4] 


1.3734 


0.4280 


2.4317 


5 


[5/5] 


1.8970 


0.6187 


2.5146 




5 


[5/5] 


1.2631 


0.4558 


2.5229 


5 


[5/6] 


1.9309 


*0.6178 


2.5105 




5 


[5/6] 


1.2435 


*0.4313 


2.4983 


5 


[6/4] 


1.9061 


0.6187 


2.6074 




5 


[6/4] 


1.2642 


0.4492 


2.7025 


5 


[6/5] 


1.9309 


*0.6267 


2.5030 




5 


[6/5] 


1.2325 


*0.4503 


2.4846 


6 


[4/6] 


1.8906 


0.6384 


2.5005 




6 


[4/6] 


1.2467 


0.5049 


2.5046 


6 


[5/5] 


1.8878 


0.6557 


2.4993 




6 


[5/5] 


1.2236 


0.6468 


2.5022 


6 


[5/6] 




*0.6091 


2.5076 




6 


[5/6] 


1.2332 


*0.3972 


2.5029 


6 


[5/7] 


*1.9425 


0.6720 


2.4446 




6 


[5/7] 


1.1881 


0.5539 


2.4788 


6 


[6/4] 


1.8897 


0.6311 


2.5098 




6 


[6/4] 


1.2428 


0.4712 


2.5305 


6 


[6/5] 


1.9346 


*0.6175 


2.5075 




6 


[6/5] 


1.2510 


*0.4210 


2.5020 


6 


[6/6] 


* 1.9248 


0.6886 


2.4268 




6 


[6/6] 


1.1799 


0.5975 


2.4780 


6 


[7/5] 


*3.0282 


0.6587 


2.4514 




6 


[7/5] 


1.1955 


0.5159 


2.4782 



Table 12: Results for the structures FFF, FFL and FFH where for the gauge parameter 
the values ^ = +2 (left table) and = +5 (right table) have been adopted. 



comparison. Actually the corresponding part of the function —g{z) is plotted in the 
vicinity oi z = 1. The ten curves in each plot correspond to ^ = —5, —2, 0, +2 and ^ = +5 
where for each value of ^ two different Fade results are shown. 

It can be seen that as for the order and a1 cases also here the choice ^ = — 2 (solid 
curve) exhibits the smoothest behaviour for z = 1 which was already expected from the 



stable behaviour of the numbers in Tabs. O and O. Thus it can be expected that for 



this values our procedure works best. However, we want to stress that — in analogy to 
the lower order analysis — the value ^ = looks very promising and also for this choice 
reasonable results can be expected. 

Adopting the choice ^ = —2 we extract from Tabs. |TT] and the following results: 

^^^^ = 2.1(4), ^^^^ = 0.8(2), = 2.4(1). (28) 

Note that the colour structure FFL demonstrates a rather strong dependence on the 
gauge parameter. To take this into account we assigned somewhat extended errors to this 
quantity in Eq. (p8|). 

Note that the results of ( ^7|) and ( ^8l) are valid for //^ = M^. With the help of the 
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n 


P.A. 


e = -5 


e--2 


C = +2 


C-+5 


4 


[1/3] 


2.4174 


2.4145 


2.4135 


2.4161 


4 


[2/2] 


2.4182 


2.4146 


2.4135 


2.4162 


4 


[3/1] 


2.4165 


2.4102 


2.4032 


2.3985 


4 


[4/0] 


2.3648 


2.3578 


2.3484 


2.3414 


5 


[1/4] 


2.4336 


2.4317 


2.4306 


2.4304 


5 


[2/3] 


2.4271 


*2.4204 


^2.4029 




5 


[3/2] 


2.4436 


2.4386 


2.4345 


2.4319 


5 


[4/1] 


2.4298 


2.4259 


2.4215 


2.4184 


5 


[5/0] 


2.3847 


2.3798 


2.3733 


2.3684 


6 


[1/5] 


2.4423 


2.4413 


2.4408 


2.4406 


6 


[2/4] 


2.4478 


2.4474 


2.4478 


2.4525 


6 


[3/3] 


2.4525 


2.4497 


2.4485 


2.4526 


6 


[4/2] 


2.4523 


2.4490 


2.4460 


2.4441 


6 


[5/1] 


2.4394 


2.4368 


2.4337 


2.4315 


6 


[6/0] 


2.3992 


2.3957 


2.3909 


2.3873 



Table 13: Results for FFH where only the low-energy moments are taken into account. 
Different choices for the gauge parameter are adopted. 



two-loop /3 1^ and three-loop 7™ 0| function it is possible to restore the ln(/i^/M^) 
terms at order a^. Thus Eq. (P reads 
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Figure 6: z dependence of —g{z) for the structures FFF, FFL and FFH for different 
values of ^. The sohd curves belong to ^ = —2. The dashed curves correspond to 
^ = —5, 0, +2, +5 (from short to long dashes). 



+ CfT^ 



1/3 _l^/2 



197 2 



(29) 



with l^M = In/i^/M^. n/ = + 1 is the number of active quark flavours. Very often the 
scale-invariant mass, /im defined through fim = "^(/^m) is of interest. Thus iterating (|29|) 
leads to 
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Inverting Eq. 
M j_ 



■m{jj,) 
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(31) 



Inn'^/m'^ln). 



If one is interested in the Eqs. (|29D , ( pO|) or ( plD the results of (|27|) and (pSD can be 
used. However, for most practical applications only the values Nc = 3, i.e. Cf = 4/3 and 
Ca = 3, and T = 1/2 are of interest. Simply adding the results of the individual colour 
factors would lead to a big overestimation of the error. A closer look to the numbers in 
Eqs. ( pTl) and (|28|) also indicates that once the numerical values for the colour factors and 
the number of light quarks are inserted numerical cancellations might occur which result 
in a loss of accuracy. Thus it is much more promising to add in a first step the results 
for the moments, of course, taking into account the proper colour factor. Afterwards the 
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n 


P.A. 




m = 1 


ni = 2 


ni =3 


n/ = 4 


n/ = 5 


5 


[4/5] 


-200.2787 


-173.6663 


-148.3787 


-124.4156 


-101.7771 


-80.4628 


5 


[4/6] 


-201.6419 


-174.8844 


-149.4553 


-125.3553 


-102.5862 


-81.1625 


5 


[5/4] 


-203.9394 


-176.7290 


-150.8970 


-126.4411 


-103.3591 


-81.6482 


5 


[5/5] 


-201.4721 


-174.7445 


-149.3422 


-125.2644 


-102.5104 


-81.0786 


5 


[5/6] 


-198.7799 


-172.8884 


-148.2001 




* - 102.2739 


* - 81.0336 


5 


[6/4] 


-202.8435 


-175.8651 


-150.2298 


-125.9387 


-102.9929 


-81.3939 


6 


[4/6] 


-201.0906 


-174.3880 


-149.0165 


-124.9749 


-102.2619 


-80.8758 


6 


[5/5] 


-200.9265 


-174.2458 


-148.8927 


-124.8668 


-102.1673 


-80.7929 


6 


[5/6] 


* - 200.4927 


* - 173.9600 


* - 148.7433 


* - 124.8358 


* - 102.2290 


* - 80.9131 


6 


[5/7] 


-200.3603 


-173.7018 


-148.3764 


-124.3940 


-101.7753 


-80.5533 


6 


[6/4] 


-201.6970 


-174.9293 


-149.4861 


-125.3673 


-102.5725 


-81.1016 


6 


[6/6] 


-200.3195 


-173.6857 


-148.3751 


-124.3879 


-101.7244 


-80.3848 


6 


[7/5] 


-202.1300 


-175.2569 


-149.7173 


-125.5125 


-102.6443 


-81.1143 



Table 14: Pade results for the sum of those contributions which don't have a closed 
heavy fermion loop, ni has been varied from to 5. 



n 


P.A. 




ni = 1 


ni = 2 


n; 3 


ni = 4 


n; = 5 


4 


[1/3] 


-0.9345 


-0.9572 


-0.9798 


-1.0024 


-1.0249 


-1.0475 


4 


[2/2] 


-0.9321 


-0.9546 


-0.9770 


-0.9995 


-1.0218 


-1.0442 


4 


[3/1] 


-0.9324 


-0.9551 


-0.9777 


-1.0003 


-1.0229 


-1.0455 


4 


[4/0] 


-0.9604 


-0.9828 


-1.0053 


-1.0277 


-1.0501 


-1.0725 


5 


[1/4] 


-0.9271 


-0.9495 


-0.9720 


-0.9944 


-1.0169 


-1.0393 


5 


[2/3] 


-0.9219 


-0.9440 


-0.9661 


-0.9882 


-1.0103 


-1.0324 


5 


[3/2] 


-0.9086 


-0.9347 


-0.9591 


-0.9827 


-1.0060 


-1.0290 


5 


[4/1] 


-0.9254 


-0.9478 


-0.9703 


-0.9927 


-1.0151 


-1.0375 


5 


[5/0] 


-0.9495 


-0.9719 


-0.9942 


-1.0166 


-1.0389 


-1.0613 


6 


[1/5] 


-0.9217 


-0.9441 


-0.9665 


-0.9888 


-1.0112 


-1.0336 


6 


[2/4] 


-0.9140 


-0.9364 


-0.9589 


-0.9813 


-1.0037 


-1.0261 


6 


[3/3] 


-0.9125 


-0.9352 


-0.9578 


-0.9803 


-1.0028 


-1.0252 


6 


[4/2] 


-0.9126 


-0.9352 


-0.9578 


-0.9803 


-1.0028 


-1.0253 


6 


[5/1] 


-0.9202 


-0.9425 


-0.9649 


-0.9872 


-1.0096 


-1.0319 


6 


[6/0] 


-0.9416 


-0.9639 


-0.9862 


-1.0085 


-1.0308 


-1.0532 



Table 15: Pade approximations performed in the variable z. No high-energy results have 
been used. Again ni has been varied from to 5, the dependence, however, is very weak. 



Pade approximations are computed. As the contributions with and without heavy quark 
loops are treated differently we consider the sum of both sets separately as a function 
of rii. In Tab. |1^ the results for the sum of the structures FFF, FFA, FFL, FAA, 
FAL and FLL is shown where rii is varied from to 5. The analog combination of the 
remaining contributions [FFH, FAH, FLH, FHH) can be found in Tab. ITB|. 



Using the results of Tabs. M and 15 with = and n^-dependent parts from Eqs. (E 
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m{M)/M 




/M 


M/'m{m) 


ni 


0{ai) 




0{ai) 


0{ai) 


0{ai) 


0(ai) 





-14.33 


-202(5) 


-11.67 


-170(5) 


13.44 


194(5) 


1 


-13.29 


-176(4) 


-10.62 


-146(4) 


12.40 


168(4) 


2 


-12.25 


-150(3) 


-9.58 


-123(3) 


11.36 


143(3) 


3 


-11.21 


-126(3) 


-8.54 


-101(3) 


10.32 


119(3) 


4 


-10.17 


-103(2) 


-7.50 


-81(2) 


9.28 


96(2) 


5 


-9.13 


-82(2) 


-6.46 


-62(2) 


8.24 


75(2) 



Table 16: Dependence of and 2;^' on rii. The choice jj? = M^, respectively, fx^ 
has been adopted. 



m 



and (|2|) the formulae (|2 

m(M) 
M 

M 

m{m) 



(|30| ) and (^) can be written in the form 
= 1 - 1.333aM + 4/ [-14.33 + 1.041n/] 



+ 4i 



-202(5) + 27.3(7)n; - 0.653n; 



1 — 1.333aM + a 



-11.67+ 1.041n/ 



+ 4/ 



-170(5) + 24.8(7)n/ - 0.653n; 



1 + 1.333a„ + [13.44 - 1.041n/] 
194(5) - 27.0(7)n/ + 0.653^ 



(32) 
(33) 

(34) 



with Om = tts (m) / vr. For simplicity /i = M and fx = m has been chosen in (^) and 
respectively. In the above equations the exact values of the nf term is displayed. Our 
method leads to 0.66(2) which is in very good agreement. This is a further justification 
for our approach. 

For completeness we also want to list the results for the values rii = 1, . . . , 5. They are 
also obtained from Tabs. |14| and ^ and summarized in Tab. |l^ where also the two-loop 
coefficients are listed. Actually the coefficients of the terms linear in ni of Eqs. (|3^)- 
(0) have been obtained by performing a fit to the three-loop results of Tab. |T6|. The 
errors of about 2-3% for the three-loop results of (p^)-(|3^) and Tab. |T6| have again been 
obtained by doubling the spread of Fade approximants. On one side this is justified with 
the behaviour at Ola"^) where the results of the first column in Tab. ^ are reproduced 
with the same order of accuracy (cf. Tab. On the other side the Fade approximants 
demonstrate more stability in the case where the moments are added and ni is fixed 
afterwards than in the case where the Fade procedure is applied to the individual colour 
structures separately. We want to stress that the errors assigned to individual terms in 
(p^)-(|3^) are, in fact, correlated. Thus for practical applications with fixed values for ni 
Tab. |T6| should be used. 
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ni 


this work 


[52 


(FAC) 


1 


52| (PMS) 


151 


1 ("large-/3o") 


2 
3 
4 
5 


143(3) 
119(3) 
96(2) 
75(2) 


152.71 
124.10 
97.729 
73.616 


153.76 
124.89 
98.259 
73.903 


137.23 
118.95 
101.98 
86.318 



Table 17: Comparison of the results obtained in this paper with estimates based of FAC, 
PMS and the "large- /3o" approximation for M/m{m). 



In our calculation we have neglected the effects due to the light quark masses. At 
order these can be taken into account by replacing the ni term in (P^), (P3|) and ( |5^ ) 
with the following expression |^ 

4 „ ^ /MA , , 

-1.041„, + - j: A ^ . (35) 

where i runs over all light quark flavours. If < r < 1 then the function A(r) may be 
conveniently approximated as follows 

A(r) = — r- 0.597 r2 + 0.230 r^ (36) 
8 

which is accurate to 1%. 



7 Discussion and Applications 

In this section we compare our results with various predictions which already exist in 
the literature. Furthermore a few applications are discussed where the relation between 
the MS and on-shell quark mass is explicitly needed. In particular we will demonstrate 
that the new term of 0{al) is indeed necessary to perform consistent analysises in some 
practical applications. 

In a first step we want to confront the estimations for the 0{a^) terms obtained with 
the help of different optimization procedures with the exact result. In ||^ the fastest 
apparent convergence (FAC) and the principle of minimal sensitivity (PMS) |Q have 
been used in order to predict the three-loop coefficient of M/m{m). In Tab. [T^the results 
are compared with ours. For rii = 2 the discrepancy amounts to only 7%. It even reduces 
to only 2% for rii = 5, i.e. in the case of the top quark. 

In Eq. (0) it is also possible to use the quantity as{M) as an expansion parameter 
(instead of as(m)). For ni = 5 the three-loop coefficient then reads 81.45(72) whereas 
the coefficient at O^a"^) does not change. For this case the authors of have obtained 



79.838 which deviates by roughly 2%. This coefficient has also been estimated in to 
be 76.17 which agrees with our central value within approximately 6%. 
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Let us compare our results also with the ones obtained in the large /3o-limit pT|, where 



/3q is the first coefficient of the QCD (3 function. The results obtained in can also 



be found in Tab. [T^. Excellent agreement below 1% is found for ni = 3. It amounts to 



roughly 5% for rii = A and 14% for rii = 5. 

In the remaining part of this section we want to discuss the application of our results 
in the context of quark production close to the threshold. Calculations in this context 
are connected to the pole mass of the involved quarks. To be specific let us consider 
the production of top quarks in e"'"e~ collisions. The corresponding physical observables 
expressed in terms of show in general a bad convergence behaviour. In the case of the 
total cross section, e.g., the next-to-next-to- leading order corrections partly exceed the 
next-to-leading ones. Furthermore the peak position which is the most striking feature of 
the total cross section and from which finally the mass value can be extracted depends 
very much on the number of terms one includes into the analysis. The commonly accepted 
explanation for this is that the pole mass is sensitive to long-distance effects which result 
to intrinsic uncertainties of order Aqcd 0, 0- other words, it is not possible to 
determine the pole mass from the analysis of the cross section at threshold with an 
accuracy better than Aqcd- 



Several strategies have been proposed to circumvent this problem [13, 14, 18]. They 
are based on the observation that the same kind of ambiguities also appear in the static 
quark potential, V{r). In the combination 2Mt + V{r), however, the infra-red sensitivity 
drops out. Thus a definition of a short-distance mass extracted from threshold quantities 
should be possible. The relation of the new mass parameter to the pole mass is used 
in order to re-parameterize the threshold phenomena. On the other hand a relation of 
the new quark mass to the MS mass must be established as it is commonly used for the 
parameterization of those quantities which are not related to the threshold. In order to do 
this consistently the three-loop relation between the MS and the on-shell mass is needed. 



In |T3[ the concept of the so-called potential mass, p^, has been introduced. Its 
connection to the pole mass is given by 

mt,ps{fj'f) = Mt - 6mt{fif) , (37) 

where Smt{fif) is connected to the potential in momentum space, V{q), through the 
equation 

SmtM = -If 7^n<r)- (38) 

Via this equation a subtracted potential, V{r,fif), is defined. The factorization scale fif 
has been introduced in order to extract the infra-red behaviour arising from the potential. 
In the combination ( |37D it cancels against the one of Mt leading to an enormous reduction 
of the long-distance uncertainties in rrit^ps Thus it is promising to formulate the 

threshold problems in terms of mt^ps[l^f) ^-nd V{r, Hf) instead of Mt and V{r). In the final 
result the dependence on /ij cancels. For the numerical analysis the value /ij = 20 GeV 



has been adopted in [T^ as its upper bound is roughly given by MtCpdsilA- 
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The computation of Smt{^f) can be performed up to 0{a^) using the exphcit results 
for V{q) |j58|. The result reads 



a. 



drutifif) = Cp/J-f — < 1 



33, 
12 ''f 



as [97 

" V [12 ~ 

9 , 305, 

— (4 / 

64^ 6 



111 

4 



121 



11 1 



r 
16 



li 

+ rii 



+ nf 



157 
324 



11/ 

54 '^^^ 



1/^ 
36 '^^ 



7145 
432 



2 r 



33623 9 



288 



13 493 
— Cs + -7^l^^f 



C2 



(39) 



a 



(6) 



(fi) and /^^ = In/ij/yU^. We are now in the position to establish a relation 



with as 

between the two short-distance masses rrit^psil^f) and mt{fi). Inserting Eqs. (|39D and (|3 
into (|37|) we get: 



mi,P5 (20 GeV) = (165.0 + 6.7 + 1.2 + 0.28) GeV , 



(40) 



where the different terms represent the contributions form order q;° to af. For the nu- 
merical values mt{mt) = 165.0 GeV and a^^\mt{mt)) = 0.1085 have been used. Note 
that the error of the 0(af) coefficient in the MS-on-shell mass relation is negligible. The 
comparison of Eq. (^0]) with the analogous expansion for Mt, 



(165.0 + 7.6 + 1.6 + 0.51) GeV, 



(41) 



shows that the potential mass can be more accurately related to the MS mass than Mt. 

The last term of the expansion in ( ^0]) is of the same order of magnitude as the error 
in the top quark mass determination at a next-linear-collider. Whereas in |l^ this term 
has been taken as uncertainty in the mass relation the error reduces significantly after 
the knowledge of the 0{al) term of the MS-on-shell relation. This can be deduced from 
the well-behaved expansion in Eq. (^Ol). The dominant error is now provided by the 
uncertainty in a^. 



A similar strategy has been proposed in [1^. There the so-called IS mass, M^^ , was 



defined as half the perturbative mass of a fictions toponium l^Si ground state which 
would exist if the top quark were stable. The relation between M^^ and Mt is given by 



IS 



Mt- 



MtClal 



e + e^^ 

TT 



2917, 275, 

+ TTTTrC2 + 



216 
193 

~6~ 



4 



97 

— + ll/^,t + ni 
6 



11 2^ 



as 

+ e-^ ( — 

TT 



1793 
12 



9 , 927, 
32^* + -' 



363,2 



+ nf 



77 1 2 

108 + 54^' + 9^' + 



1693 



72 
1/^ 



11 r - —C 



(42) 



with a. 



.(6) 



(/i) and l^^t = ln{fi/{CFasMt)). 



case of the potential mass. From the experiment the quantity M^^ is extracted, 
it has been shown that this is possible with an uncertainty of approximately 200 MeV. In 



The philosophy is very similar as in the 

In 
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a next step M}^ has to be related to the MS mass mtirrit). As the extraction of M}^ is 
based on a next-to-next-to-leading order formalism the 0{a^g) relation computed in this 
work is necessary. 

In practice one proceeds as follows: In a first step (^21) is inverted and then equated 
to (l3^). After setting = mt{fi) it is possible to solve the equation numerically for 



mt{mt). Care has to be taken in connection to the expansion parameter in Eq. (^21). 
A naive counting in powers of would lead to theoretically inconsistent results as was 



demonstrated in |]59[. A consistent treatment is obtained if an expansion parameter 



e is introduced where the terms of order are proportional to e" in Eq. (^4 ) and 
proportional to e" in Eq. p^). 

One finally arrives at the following relation between the MS and IS* mass 

mt{mt) = (175.00 - 7.60 - 0.97- 0.14) GeV . (43) 

where M^^ = 175 GeV and a^^^^Mz) = 0.118 has been adopted. Using the large- /5o 
results for the order term of ( plf ) the last term reads —0.23 [|18| which is off by more 
than 50% form the exact result. The conclusions which can be drawn from Eq. ( |^ ) are 
very similar to the ones stated above: the uncertainties due to unknown terms in the mass 
relations are negligible as compared to the error with which M^^ can be extracted from 
the experiment. The dominant uncertainty comes from the error in as which amounts for 



±0.003 to roughly 200 MeV p[ in Eq. 

For completeness we would like to mention that the introduction of short distance 
masses like rrit^ps or M}^ have significantly stabilized the position of the peak in the top 
quark cross section which is important in extracting the mass value. The normalization 
uncertainty, however, remains large when including higher order corrections. 

Also the bottom quark mass can be extracted from quantities related to the quark 
threshold. Recently [E4|, |60, EH| a precise value for the bottom quark mass has been 



determined in the context of QCD sum rules. For example, in ||2^ the on-shell mass was 
eliminated in favour of the IS mass in order to reduce the error. Once Ml^ is determined 



the analogous equation to (42) in combination with (0) can be used in order to get 
mb{mb). In the values M^^ = 4.71 ± 0.03 GeV and mb{mb) = 4.2 ± 0.06 have been 
obtained where the large-/3o approximation for the three-loop term in the MS-on-shell 
relation is part of the error in mb{mb)- Following the procedure described in |^ one 
arrives at 

mbimb) = (4.71 - 0.40 - 0.11 -0.03 ±0.03 ±0.04) GeV, (44) 

where the different terms correspond to different orders in the T-expansion. The first 
error is due to M^^ and the second one reflects the error in as{Mz) = 0.118 ±0.004 which 



is adopted from 03. Due to the nice convergent behaviour of (HI) the total error on 



mb{mb) only contains these two sources which finally leads to 

mb{mb) = 4.17 ± 0.05 GeV . (45) 
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accuracy 


0{ai) 




0{ai) 


Mb (GeV) 
mb{mb) (GeV) 


4.752 
4.362 


4.858 
4.292 


5.001 
4.322 



Table 18: Bottom quark mass determined with increasing accuracy. 



It is important to stress that taking into account of the newly computed 0{al) term in 
the MS-on-shell relation is crucial for the reliable estimation of the errors in (^). Indeed, 
a deviation of the real value for from the large-/3o estimation by, say, a factor of two, 
which one could not exclude a priori, would result to a systematic shift in mbirrih) of 
around 100 MeV. 

A somewhat different approach for the determination of both the charm and bottom 
mass has been followed in p6| . There the lower states in the heavy quarkonium spectrum 
were computed up to order a^. This allows the extraction of relatively accurate values for 
the pole masses of the charm and bottom quarks. The transformation to the MS mass 
has been performed with the help of the two-loop relation However, to the order 
the quarkonium spectrum was computed it is more consistent to use the 0{a^) relation 
provided in this paper. 

In fact, in pB| the pole mass of the bottom quark is given with a accuracy of order 



al, and a^. For the conversion to the MS scheme relation ( ^3]) is needed up to one, 
two and three loops, respectively. The same accuracy is required for the running of 
from Mz to M^. Our results are displayed in Tab. 0. As in |2§ the value af^M^) = 
0.114 has been adopted. The numerical values slightly differ from the ones given in 
which can be traced back to a minor inconsistent treatments in [^. Thus taking over 
the error estimates from 
5.001+°-^°^ 



26 1 the on-shell value for the bottom quark mass reads Mf, 



^-0.066 



GeV. It transforms to the following value for ■mb(mb) 



mbirrib) 



4.322 



+0.043 Q Y 
0.028 ^'^^ ■ 



(46) 



In a similar way also the results for the charm quark mass can be obtained. In addition 
to the bottom quark case the matching between four and five flavours has to be performed 
consistently. This means that n-loop running has to be accompanied with (n — l)-loop 
matching. For a detailed description we refer to The on-shell value 



1.8661°;}^° GeV , 



(47) 



leads to 



mATrir, 



= 1.377l°:if7GeV, 
where again the error estimates are taken over from . 



(4J 



Compared to ||2^ the inclusion of the 0(af ) terms leads to a shift in the central values 



of more than 100 MeV. In the case of the bottom quark this change is even larger than 
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the errors presented in which might indicate that the estimates were too optimistic. 
This demonstrates that a consistent treatment of the different orders in is absolutely 
crucial. 

We would like to mention that in this paper we don't intend to compare and judge the 
different methods used for extracting the bottom quark mass and the corresponding error 
estimate. For a comprehensive discussion of recent determinations of mb{mb) we refer to 

M. 



the third reference in 



8 Conclusions 

In this paper the three-loop term of order to the relation between the MS and on-shell 
definition of the quark mass is computed. This is achieved by considering expansions of 
the quark self energy for small and large external momentum followed by a conformal 
mapping and a Fade approximation. 

The three-loop computation in the limit of small momentum requires the knowledge 
of integrals which before have never been used in practical applications. We have checked 
the results available in the literature and added analytical results for the missing integrals. 

The procedure used for the computation is tested in detail at one- and two-loop order 
and compared to the known result. This also provides an estimate on the uncertainty 
which is between 2 and 3% for the order contribution. 

We discuss several important applications of the newly available term. In particular 
the relations between the MS mass and other short- distance masses are considered. It is 
shown that the three-loop term is necessary to reduce the error with which finally the MS 
quark mass can be determined. 

In conclusion we would also like to stress that our procedure is not limited to the 
determination of the conversion factor between the MS and on-shell quark masses. As a 
spin-off we also got a wealth of information about the (perturbative) quark propagator in 
QCD at order a^. It should be possible to use this approach in order to find an accurate 
semianalytic description of the propagator in all possible regions of momentum transfer 
and for a general gauge fixing condition. 
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Note added 

An analytical reevaluation of the relation between MS and on-shell quark masses have 



been reported recently in ||62[. Their result is in full agreemeent with our Eqs. (p^)-(|3 



and our Tabs. 113, IT5|, |T6| and O. It reads: 



m(M) A a. f olsV r 

^ ' - 1--— +(— j (-14.3323 + 1.0414n,) 



M 3 TT 



+ ( — (-198.7068 + 26.9239r2i - 0.65269r2f 



TT 



3 



The results of for all ten independent colour structures appearing in the a\ order are 



also in full agreement^ with ours with the only exception of the colour structure FAA. 
For this case our estimation of the error bars {z^^^ = —16.4(3) versus z^"^^ = —15.85 
as found in ||6^) seems to be sligtly too optimistic (provided, of course, that the very 
calculation of gives the correct answer for this colour structure). 
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